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Abstract. Earthquakes and aftershock sequences follow several empirical 
scaling laws: One of these laws is Bath's law for the magnitude of the largest 
aftershock. In this work, Modified Form of Bath's Law and its application 
to KOERI data have been studied. Bath's law states that the differences 
in magnitudes between mainshocks and their largest detected aftershocks are 
approximately constant, independent of the magnitudes of mainshocks and it 
is about 1.2. In the modified form of Bath's law for a given mainshock we get 
the inferred largest aftershock of this mainshock by using an extrapolation of the 
Gutenberg-Richter frequency-magnitude statistics of the aftershock sequence. To 
test the applicability of this modified law, 6 large earthquakes that occurred in 
Turkey between 1950 and 2004 with magnitudes equal to or greater than 6.9 have 
been considered. These earthquakes take place on the North Anatolian Fault 
Zone. Additionally, in this study the partitioning of energy during a mainshock- 
aftershock sequence was also calculated in two different ways. It is shown that 
most of the energy is released in the mainshock. The constancy of the differences 
in magnitudes between mainshocks and their largest aftershocks is an indication 
of scale-invariant behavior of aftershock sequences. 



1. INTRODUCTION 

An earthquake is a sudden and sometimes catastrophic movement of a part of the 
Earth's surface 1 . It is caused by the release of stress accumulated along geologic 
faults or by volcanic activity, hence the earthquakes are the Earth's natural means of 
releasing stress. When the Earth's plates move against each other, stress is put on the 
lithosphere. When this stress is strong enough, the lithosphere breaks or shifts. As 
the plates move they put forces on themselves and each other. When the force is large 
enough, the crust is forced to break. When the break occurs, the stress is released as 
energy which moves through the Earth in the form of waves, which we feel and call 
an earthquake. 

Aftershocks are earthquakes in the same region of the mainshock 1 j. Smaller 
earthquakes often occur in the immediate area of the main earthquake until the 
entire surface has reached equilibrium of stress. There are several scaling laws that 
describe the statistical properties of aftershock sequences El 0] . Gutenberg-Richter 
frequency-magnitude scaling law is widely known by seismologists and scientists. On 
the Richter scale, the magnitude (M) of an earthquake is proportional to the log of 
the maximum amplitude of the earth's motion. What this means is that if the Earth 
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moves one millimeter in a magnitude 2.0 earthquake, it will move ten millimeters in 
a magnitude 3.0 earthquake, 100 millimeters in a magnitude 4.0 earthquake and ten 
meters in a magnitude 6.0 earthquake. So, the amplitude of the waves increases by 
powers of ten in relation to the Richter magnitude numbers. Therefore, if we hear 
about a magnitude 8.0 earthquake and a magnitude 4 earthquake, we know that the 
ground is moving 10, 000 times more in the magnitude 8.0 earthquake than in the 
magnitude 4.0 earthquake. Numbers for the Richter scale range from to 9.0, though 
no real upper limit exists. The difference in energies is even greater. For each factor 
of ten in amplitude, the energy grows by a factor of 32. When seismologists started 
measuring the magnitudes of earthquakes, they found that there were a lot more small 
earthquakes than large ones. Seismologists have found that the number of earthquakes 
is proportional to 10~ M/ . They call this law "The Gutenberg-Richter Law" [HI EI EJ- 
In seismological studies, the Omori law, proposed by Omori in 1894, is one of the 
few basic empirical laws 6 . This law describes the decay of aftershock activity with 
time. Omori Law and its modified forms have been used widely as a fundamental 
tool for studying aftershocks 0. Omori published his work on the aftershocks of 
earthquakes, in which he stated that aftershock frequency decreases by roughly the 
reciprocal of time after the main shock. An extension of the modified Omori's law is 
the epidemic type of aftershock sequences (ETAS) model 0- It is a stochastic version 
of the modified Omori law. In the ETAS model, the rate of aftershock occurrence is 
an effect of combined rates of all secondary aftershock subsequences produced by each 
aftershock jHUH]- 

The third scaling law relating the aftershocks is Bath's law. The empirical Bath's 
law states that the difference in magnitude between a mainshock and its largest 
aftershock is constant, regardless of the mainshock magnitude and it is about 1.2 

Ei m h nni - That is 

Am = m ms - m2 ax (1) 

with m ms the magnitude of the mainshock, m™ s ax the magnitude of the largest 
detected aftershock, and Am approximately a constant. 

In this article we study the modified form of Bath's law 0). To study 
the aftershock sequence in the North Anatolian Fault Zone (NAFZ) we get the 
largest aftershock from an extrapolation of the G-R frequency-magnitude scaling 
of all measured aftershocks. We test the applicability of Bath's law for 6 large 
earthquakes on the North Anatolian Fault Zone (NAFZ). The emprical form of 
Bath's law states that the difference magnitude between a mainshock and its largest 
aftershock is constant, independent of the magnitudes of mainshocks. We also analyze 
the partitioning of energy during a mainshock-aftershock sequence and its relation to 
the modified Bath's law. 

2. BATH'S LAW AND ITS MODIFIED FORM 

Bath's law states that the differences in magnitudes between mainshocks and their 
largest aftershocks are approximately constant, independent of the magnitudes of 
mainshocks. In modified form of Bath's law for a given mainshock we get the 
inferred largest aftershock of this mainshock by using an extrapolation of the 
Gutenberg-Richter frequency-magnitude statistics of the aftershock sequence. The size 
distribution of earthquakes has been found to show a power law behavior. Gutenberg 
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and Richter, introduced the common description of the frequency of earthquakes: |S] 

logigN(> m) = a — bm (2) 

where iV(> m) is the cumulative number of earthquakes with magnitudes greater 
than m occurring in a specified area and time window. On this equation a and b are 
constants. This relation is valid for earthquakes with magnitudes above some lower 
cutoff m c . Earlier studies |3J0J^3 gave an estimate for this "b" value between 0.8 and 
1.2. The constant "a" shows the regional level of seismicity and gives the logarithm of 
the number of earthquakes with magnitudes greater than zero [13 E]- In our analysis 
a-value is in the range 3.8 < a < 6.5. Aftershocks related with a mainshock also 
satisfy G-R scaling (j2J to a good approximation 3 . In this case N(> to) is the 
cumulative number of aftershocks of a given mainshock with magnitudes greater than 
m. We offer to extrapolate G-R scaling (01 for aftershocks. Our aim is to obtain an 
upper cutoff magnitude in a given aftershock sequence. We find the magnitude of 
this inferred "largest" aftershock m* by formally taking N(> to*) = 1 for a given 
aftershock sequence. 

a = bm* (3) 

This extrapolated m* value will have a mean value and a standard deviation from the 
mean value. We apply the Bath's law to the inferred values of m* and then, we can 
write 

Am* — m ms — m* (4) 

where m ms is the magnitude of the mainshock and Am* is approximately a constant. 
Substitution of equations J2J and (@J into equation J2J) gives 

logw[N(> to)] = b(m ms — Am* — to) (5) 

with b, m msi and Am* specified, the frequency-magnitude distribution of aftershocks 
can be determined using equation In extrapolating the G-R scaling (j2J the slope 
of this scaling or b-value plays an important role in estimating the largest inferred 
magnitude to*. 

3. APPLICATION OF THE MODIFIED FORM OF BATH'S LAW TO 
THE NORTH ANATOLIAN FAULT ZONE (NAFZ) 

We applied modified form of Bath's law by considering 6 large earthquakes on the 
NAFZ. These earthquakes occurred between 1950 and 2004. The data were provided 
by Bogazici University Kandilli Observatory and Earthquake Research Institute 
(KOERI) 12 . The 6 earthquakes considered had magnitudes m ms > 6.9. The 
important point is that they were sufficiently separated in space and time so that 
no aftershock sequences overlapped with other mainshocks. Earthquakes form a 
hierarchical structure in space and time. Therefore, in some cases it is possible to 
discriminate foreshocks, mainshocks, and aftershocks. But, generally this classification 
is not well defined and can be ambiguous. One of our main problems in the study 
of aftershocks is to identify what is and what is not an aftershock JSj- To specify 
aftershocks we defined space and time windows for each sequence. In each case we 
consider a square area centered on the mainshock epicenter. The linear size of the box 
is taken to be of the order of the linear extent of the aftershock zone L, which scales 
with the magnitude of the mainshock as 

L = 0.02 x 10 a5mms A:TO (6) 
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Figure 1. Frequency-magnitude distribution of Kocaeli-Golciik earthquake 



This equation was given by Yan Y. Kagan ^1]. Previously, Shcherbakov and Turcotte 
used the same scaling arguments for 10 earthquakes in California [HIE]- Time intervals 
of 92, 183, 365, 730, and 1095 days are taken except Canakkale-Yenice, Mu§-Varto, and 
Adapazari-Mudurnu earthquakes. It should also be noted that for all 6 earthquakes 
we took mr,, Richter magnitudes. 

For Kocaeli-Golciik earthquake the mainshock is m ms = 7.4 and the largest 
detected aftershock had a magnitude m™ aa: = 5.8. From equation the difference 
in magnitude between the mainshock and largest aftershock is Ato = 1.6. We have 
correlated the aftershock frequency magnitude data given in Figure 1 with G-R scaling 
© and find b = 0.91 ± 0.05 and a = 4.97 ± 0.20. From equation J3J the inferred 
magnitude of the largest aftershock is to* = 5.46 ± 0.37. From equation Q the 
difference in magnitude between the mainshock and the inferred largest aftershock is 
Ato* = 1.94±0.37. We applied the same procedure to the all 6 earthquakes and found 
a, b, to*, and Am* parameters. 

For these 6 earthquakes, the a, b, m ms , m™ aa: , Am, to* , and Am* values are 
given in Tabled 

According to our results, the mean of the differences between mainshock and 
largest detected aftershock magnitudes is Ato = 1.63 ±0.23. The mean of the inferred 
values of Ato* obtained from the best fit of equation (0 is Ato* = 1.42 ± 0.18. In 
addition for these earthquakes the mean of b values is b = 0.81 ± 0.02. 
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Table 1. Summaries of the data and results that show b and a parameters 



Earthquake 


Date (mm/dd/yy) 


b 


a 


Canakkalc-Ycnicc 


03/18/53 


0.82 ±0.07 


4.69 ±0.33 


Bolu-Abant 


05/26/57 


0.61 ±0.05 


3.81 ±0.24 


Mu§-Varto 


08/19/66 


0.95 ±0.09 


5.20 ±0.37 


Adapazari-Mudurnu 


07/22/67 


1.22 ±0.10 


6.62 ± 0.45 


Kocacli-Golciik 


08/17/99 


0.91 ±0.05 


4.97 ± 0.20 


Diizce 


11/12/99 


0.80 ±0.03 


4.85 ±0.13 



Table 2. Summaries of the data and results that show m ms , m™ s ax , Am, m*, 
Am* parameters 



Earthquake 




nb as 


Am 


m* 


Am* 


Canakkale-Yenice 


7.2 


5.4 


1.8 


5.72 ±0.63 


1.48 ±0.63 


Bolu-Abant 


7.1 


5.9 


1.2 


6.25 ±0.65 


0.85 ±0.65 


Mu§-Varto 


6.9 


5.3 


1.6 


5.47 ±0.65 


1.43 ±0.65 


Adapazari-Mudurnu 


7.2 


5.4 


1.8 


5.43 ±0.58 


1.77 ±0.58 


Kocacli-Golciik 


7.4 


5.8 


1.6 


5.46 ±0.37 


1.94 ±0.37 


Diizce 


7.2 


5.4 


1.8 


6.06 ±0.28 


1.14 ±0.28 
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Figure 2. Dispersion of the magnitude differences Am and Am* on the 
mainshock magnitude m ms . White symbols correspond to Am values and black 
symbols correspond to Am* values. 
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Figure 3. The relation between Am and Am*. Line 1 shows the harmony of 
our data with the Bath's Law. Line 2 corresponds to y = x line and shows the 
harmony of our data with the Modified Form of Bath's Law. 



4. Radiated Energy During an Earthquake 

Seismologists have more recently developed a standard magnitude scale that is called 
the moment magnitude, and it comes from the seismic moment. To understand the 
seismic moment, we need to go back to the definition of torque. A torque is an agent 
that changes the angular momentum of a system. It is defined as the force times the 
distance from the center of rotation. Earthquakes are caused by internal torques, from 
the interactions of different blocks of the earth on opposite sides of faults. It can be 
shown that the moment of an earthquake is simply expressed by 

M = iiAd (7) 

where M =Moment, /z=Rock Rigidity, A=Fault Area, and d=Slip Distance. Both the 
magnitude and the seismic moment are related to the amount of energy that is radiated 
by an earthquake. Radiated energy is a particularly important aspect of earthquake 
behavior, because it causes all the damage and loss of life, and additionally, it is the 
greatest source of observational data. So, the seismic radiated energy is an important 
physical parameter to study on earthquakes. The relationships between the radiated 
energy, stress drop, and earthquake size provides information about the physics of the 
rupture process. Richter and Gutenberg, developed a relationship between magnitude 
and energy. Their relationship is: 

log w [E(m)] = |m + 11.8 (8) 

It should be noted that in this relation E(m) is not the total "intrinsic" energy of the 
earthquake. It is only the radiated energy from the earthquake and a small fraction of 
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the total energy transferred during the earthquake process. We can write this equation 
in this form ■ 

3 

log w [E(m)} = -m + log w E (9) 

with Eq = 6.3 x 10 J. Our aim is to determine the ratio of the total seismic energy 
radiated in the aftershock sequence to the seismic energy radiated in the mainshock. 
This relation can be used directly to relate the radiated energy from the mainshock 
E rns to the moment magnitude of the mainshock m ms 

E ms = E 10 {3/2)m - (10) 

4--1- The First Calculation Method To Find The Energy Ratio Between Mainshock 
and Aftershock Sequences 

The total radiated energy in the aftershock sequence E as is obtained by integrating 
over the distributions of aftershocks pjj. This can be written 

E as = / E (m)(- — )dm (11) 

Taking the derivative of equation J2J with respect to the aftershock magnitude m we 
have 

dN = -b(lnl0)10 a - bm dm (12) 
Putting equation (|12fl into equation (|11|) gives 

/ m T7iaa; 
E(m)10- bm dm (13) 
-OO 

In addition, if we turn back to equation {TUJ) and put it to equation l|13f) we get 

lQ(3/2-b)m dm ( U ) 

-OO 

Then we take this integral and we find 

E as = (3 ! b 2b) £olO a 10 (3/2 ~ fc)mrX (15) 

To find the ratio of the total radiated energy in aftershocks E as to the radiated energy 
in the mainshock E ms , we divide equation l|15H to equation IjlOl) . Then we get the 
result 

Jel = 26 lono-fr^io- 3 / 2 ^-™^*) (16) 

E ms (3 - 26) 

We know that Am — m ms — m™ s ax so equation l|16|) takes this form: 

— = - 26 iono- b <° x io- 3 / 2 ( Am ) (17) 

E ms (3 - 26) 

From the equation (|17|D. the fraction of the total energy associated with aftershocks is 
given by 

Eas _ 1 f 18) 

Ems + E as 1+ (3_^)10 3 / 2Am 10-( a - bm - *) 1 ^ 

For the 6 earthquakes considered in the previous section we had put the b, a, Am and 
m™" values to equation (|TH|) individually. Our aim is to find E E ^ E values for 6 
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Table 3. Summaries of the data and results for energy values that were taken 
from the first energy calculation method 



Earthquake 


a 


b 


ln as 


Am 


E m „+E„.„ 


Canakkale-Yenice 


4.69 ±0.33 


0.82 ±0.07 


5.4 


1.8 


0.004 


Bolu-Abant 


3.81 ±0.24 


0.61 ±0.05 


5.9 


1.2 


0.017 


Mu§-Varto 


5.20 ±0.37 


0.95 ±0.09 


5.3 


1.6 


0.010 


Adapazan-Mudurnu 


6.62 ±0.45 


1.22 ±0.10 


5.4 


1.8 


0.009 


Kocaeli-Golcuk 


4.97 ±0.20 


0.91 ±0.05 


5.8 


1.6 


0.002 


Diizce 


4.85 ±0.13 


0.80 ±0.03 


5.4 


1.8 


0.008 



earthquakes. 

According to our results, we find the mean energy E E ^ a E — = 0.008 with a 
standard deviation er-g = 0.005. Consequently, we find that for these earthquakes 
on average about 99.2 per cent of the available elastic energy is released during the 
mainshock and about 0.8 per cent of energy is released during the aftershocks. 

4-2. The Second Calculation Method To Find The Energy Ratio Between Mainshock 
and Aftershock Sequences 

Additionally, from the study of Shcherbakov and Turcotte in 2004 HUE], we may derive 
the same energy ratio in terms of b and Am* values. 

The total radiated energy in the aftershock sequence E as is obtained by 
integrating over the distributions of aftershocks j2j. This can be written 



E„ 



E(m)( — - — )dm 
dm 



(19) 



Taking the derivative of equation (Jjjjwith respect to the aftershock magnitude m we 
have 



dN = -b(lnl0)10 b( - m ™- Am *- m Um 
Putting equation (|20fl into equation (|19|l gives 

/tn* 
E(m)l0- bm dm 
-OQ 

In addition, if we turn back to equation l|10[) and put it to equation l|21|l we get 

E as = b(lnl0)10 b ( m ™°- Am '^E / 10 (3/2 - b)m dm (22) 



(20) 
(21) 



Then we take this integral and we find 

2b 



E„ 



(3 - 2b) 
Using equation (0} we find 

2b 
(3 - 2b) 



-EnlOW 2-6 ' 1 ™* ^Q6(m m »-Am*) 



E ns — 



E 10 3 / 2( - r 



,-Am*) 



(23) 



(24) 
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Table 4. Summaries of the data and results for energy values that were taken 
from the second energy calculation method 



Earthquake 


b 


Am* 




Canakkale-Yenice 


0.82 ±0.07 


1.48 ±0.63 


0.007 


Bolu-Abant 


0.61 ±0.05 


0.85 ±0.65 


0.035 


Mu§-Varto 


0.95 ±0.09 


1.43 ±0.65 


0.012 


Adapazan-Mudurnu 


1.22 ±0.10 


1.77 ±0.57 


0.010 


Kocaeli-Golciik 


0.91 ±0.05 


1.94 ±0.37 


0.002 


Diizce 


0.80 ±0.03 


1.14 ±0.28 


0.022 



To find the ratio of the total radiated energy in aftershocks E as to the radiated energy 
in the mainshock E ms , we divide equation 1)24(1 to equation 1)10(1 . Then we get the 
result 



Ea^ = M 3 /2Am- (25) 

E ms (3-26) 1 j 

If we further assume that all earthquakes have the same seismic efficiency (ratio of 
radiated energy to the total drop in stored elastic energy), then this ratio is also the 
ratio of the drop in stored elastic energy due to the aftershocks to the drop in stored 
elastic energy due to the mainshock. From equation l(25|) the fraction of the total 
energy associated with aftershocks is given by 



2b 



E as 



1 



E n 



E„ 



1 



2b 



1Q 3/2Ar 



(26) 



For the 6 earthquakes considered in the previous section we had put the b and 
m* values to equation ((261) individually. Our aim is to find E E ^ E — values for 6 
earthquakes considered. The obtained results are summarized in Table Q. 

According to our results, we find the mean energy E E ^ E — = 0.015 with a 
standard deviation er-g = 0.012. Consequently, we find that the ratio of radiated 
energy in aftershocks to the radiated energy in the mainshock is constant. This is 
consistent with the generally accepted condition of self-similarity for earthquakes. For 
these earthquakes on average about 98.5 per cent of the available elastic energy goes 
into the mainshock and about 1.5 per cent into the aftershocks. 



5. Discussion 



Earthquakes occur in clusters. After one earthquake happens, we usually see others 
at nearby or identical location. Clustering of earthquakes usually occurs near the 
location of the mainshock. The stress on the mainshock's fault changes drastically 
during the mainshock and that fault produces most of the aftershocks. This causes a 
change in the regional stress, the size of which decreases rapidly with distance from 
the mainshock. Sometimes the change in stress caused by the mainshock is great 
enough to trigger aftershocks on other, nearby faults. It is accepted that aftershocks 
are caused by stress transfer during an earthquake. When an earthquake occurs there 
are adjacent regions where the stress is increased. The relaxation of these stresses 
causes aftershocks El El El El H3 ED 
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Table 5. Comparision the results 



Parameters 


Turcotte and Shcherbakov 


Kurnaz and Yalcin 


Am 


1.16 ±0.46 


1.63 ±0.23 


Am* 


1.11 ±0.29 


1.42 ±0.18 


Has 

E m „+E„„ 


0.038 


0.015 



Several scaling laws are also found to be universally valid for aftershocks El EI] • 
These are: 

(i) Gutenberg-Richter frequency-magnitude scaling 

(ii) The modified Omori's law for the temporal decay of aftershocks 

(iii) Bath's law for the magnitude of the largest aftershock 

In this work we used both Bath's law and G-R scaling. Our aim is to find an upper 
cutoff magnitude m* for a given aftershock sequence. Using relation @, we get related 
a and b values in the G-R scaling. Bath's law states that, to a good approximation, 
the difference in magnitude between mainshock and its largest aftershock is a constant 
independent of the mainshock magnitude. A modified form of Bath's law was proposed 
by Shcherbakov and Turcotte in 2004 00]. They considered 10 large earthquakes 
that occurred in California between 1987 and 2003 with magnitudes equal to or 
greater than m ms > 5.5. According to their theory the mean difference in magnitudes 
between these mainshocks and their largest detected aftershocks is 1.16 ± 0.46. This 
result is consistent with Bath's Law. They found the mean difference in magnitudes 
between the mainshocks and their largest inferred aftershocks is 1.11 ± 0.29. They 
also calculated the partitioning of energy during a mainshock-aftershock sequence and 
found that about 96 per cent of the energy dissipated in a sequence is associated with 
the mainshock and the rest (4 per cent) is due to aftershocks. Their results are given 
in Table El We applied the Modified Form of Bath's Law to our 6 large earthquakes 
that occurred on the North Anatolian Fault Zone (NAFZ) in Turkey. We followed the 
same calculation process. 

According to TableEl for the North Anatolian Fault Zone (NAFZ), a large fraction 
of the accumulated energy is released in the mainshock and only a relatively small 
fraction of the accumulated energy is released in the aftershock sequence. The results 
of Turcotte and Shcherbakov are for the ten earthquakes in California on the San 
Andreas Fault Zone. Although SAFZ (in California) and NAFZ (in Turkey) have the 
same seismic properties, the released energy during the mainshocks in the NAFZ is 
much greater than the released energy during the mainshocks in the SAFZ. 

Figure 7 shows the dispersion of the magnitude differences Am and Am* on the 
mainshock magnitude m ms . In this figure, white symbols correspond to Am values 
and black symbols correspond to Am* values. Square, circle and triangle were used 
to prevent the coincides of the data on the figure; because for six earthquakes, three 
of them have the same mainshock magnitude m ms = 7.2. 

Additionally, Figure 8 gives us the relation between Am and Am*. In this figure, 
line 1 shows the harmony of our data with the Bath's Law. According to Figure 8, 
our data do not show harmony with the Bath's Law. Bath's Law states that the 
difference in magnitude between a mainshock and its largest detected aftershock is 
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constant, regardless of the mainshock magnitude and it is about 1.2 0J [7J 11U| . 
But in Figure 8, only one earthquake has Am values equal to 1.2. The other five 
earthquakes have Am values greater than 1.2. Consequently, only 17 per cent of our 
data show harmony with the Bath's Law. The rest part ( 83 per cent) of our data do 
not show harmony with the Bath's Law. 

The constancy of the differences in magnitudes between mainshocks and their 
largest aftershocks is an indication of scale-invariant behavior of aftershock sequences. 

In Figure 8, line 2 shows the harmony of our data with the Modified Form of 
Bath's Law. This line corresponds to y = x line. If the Modified Form of Bath's Law 
gave us perfect results, Am and Am* values would be close to each other along this 
line. Hence, they would be the near of line 2. But in Figure 8, only one earthquake 
takes place on the upper side of this line. The remaining five earthquakes take place 
on the lower side of this line. Consequently, our data do not show harmony with the 
Modified Form of Bath's Law. 

The other important conclusion is that we know most of the energy is released 
during the mainshock. Therefore, after the mainshock the community and government 
may begin their work to rescue people from the debris without wasting any time. 
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